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A method is described for calculating corrections to the Boltzmann/Chapman-Enskog 
^^ I analysis of lattice gases due to the buildup of correlations. It is shown that renormalized 

transport coefficients can be calculated perturbatively by summing terms in an infinite series. 
A diagrammatic notation for the terms in this series is given, in analogy with the Feynman 
diagrams of quantum field theory. This theory is applied to an example lattice gas and 
shown to correctly predict experimental deviation from the Boltzmann prediction. 



I. INTRODUCTION 



The central problem in the theoretical analysis of a lattice gas is the determination of the macroscopic 
hydrodynamic equations obeyed by the conserved quantities of the system, with transport coefficients ex- 
pressed as functions of those conserved quantities. Most analyses of this problem make use of the Boltzmann 
molecular chaos assumption. This assumption, which neglects correlations between colliding particles, makes 
possible the derivation of a Boltzmann equation for the single-particle distribution function. For lattice gases 
k>( \ satisfying semi-detailed balance, this Boltzmann equation has Fermi-Dirac equilibria. The Chapman-Enskog 

'V^ ■ procedure of linearizing about these equilibria then yields the hydrodynamic equations satisfied by the sys- 

C^ I tem. The transport coefficients that appear in these equations are related to the eigenvalues of the linearized 

Boltzmann collision operator at equilibrium, and these are calculable as functions of the conserved quantities. 
Unfortunately, the first step in this analysis - the Boltzmann molecular chaos assumption - iaxiften invalid, 
and its use can lead to serious errors in the transport coefficients, as has been shown empiricallyEJ. To see why, 
recall that the transport coefficients measure the rate at which the propagation/collision process smoothes 
gradients. These same collisions, however, generate correlations among the bits of the lattice gas, and these 
correlations can propagate about and subsequently correct the Boltzmann collision operator. Thus, in an 
exact theory, the renormalized linearized collision operator is the sum of a bare Boltzmann collision operator, 
and a series of corrections due to the presence of correlations. 
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For example, consider an event in which two particles emerge from a collision (and thereby acquire a 
correlation), move about in a background of uncorrelated particles, and later recoUide. This process can be 
thought of as a one-loop correction to the Boltzmann approximation. Further refinements can be obtained 
by including more intricate diagrams - with multiple loops, nested loops, etc. - to account for the interaction 
of the correlated quantities with the background. Standard field-theoretic techniques can then be used to 
approximate these diagrammatic sums. 

In this work, we outline a systematic approach to the exact analysis of a general lattice gas. We show how 
to construct a diagrammatic series for the transport coefficients of any lattice gas satisfying semi-detailed 
balance^, and we provide closed-form expressions for the vertex coefhcients. These results constitute the 
central part of this paper and are contained in Section 3. The methods that we use here are analogous to 
the cluster expansion techniques that have been used for many years in the analysis of continuum fluidsEI. 
The discrete nature of lattice gases, however, necessitates a construction which differs substantially from the 
continuum analysis. 

The formalism described here is exact, in that it yields the transport coefficients with no approximation. 
Of course, it is usually a very difficult combinatorial task to sum the diagrammatic series involved, and, 
more often than not, one must content oneself with an approximation to this series. In doing so, however, it 
is possible to make contact with other methods of computing this correction. In particular, it is possible to 
understand kinetic ring theory, and BBGKY truncations in terms of truncations of the exact diagrammatic 
series, and to characterize the classes of diagrams that they retain, and those that they ignore. 

Our analysis differs from previous lattice gas analyses in a few important ways. We formulate the theory 
in terms of connected correlation functions. A similar expansion in terms of products of fluctuations was 
previously used in&a and inQu, where the ring kinetic theory was derived for lattice gases, and for lattice 
Lorentz gases, respectively. The use of the connected correlation functions simplifies the form of the complete 
diagrammatic expansion for the kinetic theory significantly. Previous work on the kinetic ring theory of lattice 
gases has used the Green-Kubo formalism to obtain the series for the transport coefficients. In this paper, 
we use the Chapman-Enskog theory instead. In this way, we get the Boltzmann approximation at zeroth 
order. All higher-order terms in our series are thus corrections to the Boltzmann approximation. In the 
Green-Kubo theory, by contrast, it is necessary to sum an infinite number of terms just to get the Boltzmann 
approximation. 

We describe in Section 4 a variety of simplifications and approximations which can be made to simplify the 
numerical computation of the renormalized transport coefficients. The simplest of these approximations is 
the ring approximation, in which correlations between more than 2 particles are neglected, and correlations 
are not allowed to interact. Thus, the results ofl are described in our theory by the restriction of the 
diagrammatic sum to the simple set of ring diagrams. 

Finally, to make all this more concrete, in Section 5 a worked example is given, and the results are 
compared with experiment. 

II. DEFINITIONS AND NOTATION 

A lattice gas is generally described by a state space and a time-development rule. The state space is 
defined by associating n bits with each point on a lattice L. (Bits are variables taking values in {0, 1}.) We 
define the set of bits at a general lattice site to be B, so that \B\ = n. We denote the total number of bits on 
the lattice by N — n\L\. For each value of the discrete time parameter t, we write the values of the bits as 
n*(x, t), where x G L and i € B. These bits can be thought of as a set of occupation numbers for individual 
particle states. Each of the bits n* is associated with a lattice vector c*, which represents the velocity vector 
of the associated particle. 

We define the set S of possible states of the bits at a general lattice point at a fixed value of t to be 



^at lowest order 



S = {s : s C B}, 

where a state s is associated with the set of bits taking the value 1. Note that \S\ = 2". 

We shall often want to refer to the N = n\L\ bits of the lattice in a uniform fashion, so we introduce 
an enumeration of these N bits, given by a 1-1 correspondence between the sets B = {1,2,...,A^} and 
B X L = {(i,x) : i G _B,x G L}. In this notation, a single bit of the lattice gas is written as n°, a € B. To 
relate this notation to the more explicit (i,x) notation, we express the above 1-1 correspondence by writing 
a and («,x) as functions of one another, so that 

n^it) =n*(")(x(a),i) 

and 

We shall use both notations interchangeably throughout this paper. 

III. THE RENORMALIZED CHAPMAN-ENSKOG ANALYSIS 
A. Ensemble averages 

We will denote the mean of an arbitrary product of the n^'s by 

7V" = mnA, aC{l,...,N}. 

\aea I 

There are 2^ — 1 independent N"^s, since N'^ = (1) = 1. We will sometimes use a roman index to denote 
an index set with a single element, as in N"' = N^'^^ 

We introduce an advection operator, A'^o, which is a permutation matrix on the 2^ — 1 dimensional space 
of means. A" a takes a set of particles (3 = {&i, . . . , bq} to the advected set of particles a = {ai, . . . , Uq}, 
where aj = a{i{bj),x{bj) + c'^^j^). Note that this generalizes the usual concept of an advection operator, 
that propagates single particles from site to site, to one that can propagate arbitrary multipoint quantities 
(e.g., connected correlation functions) from site to site. 

B. Boltzmann Analysis 

The evolution of a lattice gas in one time step can be divided into two parts. The first part is a collision 
phase in which the n bits at each point alter their values in some specified way. The second part is a 
propagation phase in which the value of each bit moves to the corresponding bit of the lattice point to which 
it is carried by its associated velocity vector. The complete time-development equation is then written as 

?i*(x-f c',t + At) = n^(x,t)+tj^(7i*(x,t)). (III.l) 

The collision operator uj^ is generally a nonlinear function of its arguments. Note that the form f{z*) is used 
to indicate that a function / depends on a quantity z for all possible values of the index *. 
The Boltzmann molecular chaos assumption can be written 

(c^^K(x,t)))«c.*((n*(x,t))). 

This is inexact when w' is nonlinear. 

Taking the ensemble average of Eq. ( III.l ), and making the molecular chaos assumption leads to the 
Boltzmann equation. 



N\t + At) = A\ (n'' (i) + n''^"^ fjsiai*Mc)) (^)A 



(III.2) 



where repeated indices (c in this case) are summed. This is a closed system of equations for the A^^'s. 

The usual Chapman-Enskog analysis proceeds by expanding this equation of motion using Navier-Stokes 
ordering: A formal expansion parameter e is introduced, such that |c*| ~ e. The time step is ordered as 
At ~ e^, as is appropriate for viscous or diffusive processes. Since the lattice gas is assumed to satisfy semi- 
detailed balance, there is a Fermi-Dirac equilibrium. It is assumed that deviations from this equilibrium are 
of order e. We will denote the equilibrium values of the single-particle means by Nq. 

Expanding the collision operator in (III. 2) to first order in e gives 



N^{t + At)^A\(N''{t) + 



eJ*^'^^ ^"(j'''(c)) 



it) 



(III.3) 



where J*- = d^^{N*)/dN^\ is the Jacobian matrix of the collision operator evaluated at the equilibrium. 

The rest of the Chapman-Enskog analysis consists of associating the hydrodynamic and kinetic modes of 
the system with eigenvectors of J, and writing the hydrodynamic equations in such a way that the transport 
coefficients are functions of the eigenvalues of J. For a more in-depth review of the details of this procedure, 
seeQ. 



C. Exact Analysis 



We will now drop the Boltzmann assumption and show that the exact equation of motion leads in the 
scaling limit to an equation identical to (III. 3), but with a renormalized J matrix. We will expand about 
the same equilibrium as in the usual analysis. For any lattice gas satisfying semi-detailed balance, this 
equilibrium exists; in the vicinity of this equilibrium, correlations of at most order e are generated in a single 
time step. 



We begin with the exact time-development equation ( 111.1). By taking the ensemble average of the product 
of this equation over all a in an arbitrary set a C B, we can write the exact equation for an arbitrary 
multipoint mean at time t + At in terms of multipoint means at time t. We have 



iv"(i+Ai)=/n^"(i+At))-5]^"^(n[n'w+^'^''K(x(6),t)) 



\bei3 



(III.4) 



To express the right-hand side in terms of multipoint means, it will be convenient to rewrite this equation 
in a more compact notation. For a set (3 Q B, let us define Lp to be the subset of points in L which contain 
at least one particle in the set f3; that is, 

Lf^ = {y E L : x(6) — y for some b £ P}. 

Similarly, we define /3x to be the set of i's corresponding to the particles in f3 at the point x; that is, 

/3x = {i e S : a{i,x) G /3}. 



We can now factorize the product appearing on the right-hand side of ([II.4) into contributions from each 
of the points in Lp, by writing 



[];p(t)+u;^W(n*(x(6),t))' 



fce/3 



n n [«^(x,t)+^^(n*(x,i))] 



The innermost product on the right now depends only on quantities at a single site, x. 
The functions lu^ (n* ) can be expressed as polynomials in the n' 's of the form 



L.'(n*)=^fcVn^ 



yCB 



jev 



where the fc*^, are coefHcients which may depend only on stochastic variables at each lattice site, and which 
are constant for deterministic lattice gases. Thus, we can write 



n[n' + c.Xn*)]=^^;^n^ 



iG/J 



vCLB jev 



where the quantities v'^j^ may contain stochastic variables at each site. Taking the ensemble average over 
any such variables, we get the mean vertex coefficients V^^, 



V 



(III.5) 



The lattice gas state transition probabilities A{s — > s') may be interpreted as elements of a collision matrix 
on the space of probabilities, P*, in the sense that the post-collision probability of a state s' is given by 



J2Ais^s')P'' 



(III.6) 



Similarly, the matrix V^^, can be interpreted as a collision matrix on the space of means. Since the proba- 
bilities and the multipoint means are straightforwardly related, it is possible to relate the matrix V to the 
transition probabilities. We get 






(III.7) 



where we have identified the state s with the set of bits which are 1 in that state (s C B). Clearly, the 
2^" matrix elements V^, depend only on the sets /i and v, and on the form of the collision operator. In 

particular, they do not depend on x, or on the values of the n^'s. Note that V — S'^ , regardless of the 
specific latti ce gas or collision rule. 

Eq. (111.4) can now be rewritten in the form 



N''{t + At)=A''i3K''\N^t), 



(HI. 

where K is an operator expressing the complete collision part of the time development, given by 

-/3 _ TT T^/3: 



K"- 



Uv' 



1 11 7x 



(III.9) 



We shall now rewrite the dynamical equation ( III. 8 ) in the language of connected correlation functions 
(CCF's). For each a C {1, . . . , TV}, there is a CCF, written F". It is simplest to relate the CCF's to the 
means by the equation 



7r(Q) 



(III.IO) 



where 7r(a) is the set of all partitions of a into disjoint subsets, Ci; ■ ■ • jC?- Thus, for example, we have 
N" = T" and 



N 



abc 



^abc 



aj^bc I T^bj^ac 



T^aj^i 



j^oj^a 



cj^ab 



pep 



T^ar^&T^c 



For notational convenience, an index set with a circumflex, e.g. d, is assumed to have at least two elements. 
We refer to the me an (CCF) of a set with q elements as a q-mean (g-CCF). 

Equation (III.IC) can be inverted by inducti on on |a| to get an expression of the form F" — g^lN*), where 
g = /^^. We can now use / and g to convert (111.8) into a dynamical equation for the CCF's, giving 



F"(t + At) = ^%$'3(r*), 



(III.ll) 



where 



we have used the fact that A is a permutation matrix, so that g and A commute. 

For an ensemble of independent variables, all CCF's except the 1-CCF's vanish. The equilibrium distri- 
bution is such an ensemble, so we can take any term F" to be of order e. In the asymptotic analysis, only 
first order terms in e need be kept in the expression <I>''(F*). Thus, the dynamical equation for the CCF's 
can be linearized to 



F"(i + At) = A"^ (eKf iVf (t) + Kf^r^it) 

P \ a 7 



where 



/C' 



a$/3 



7 dY-y 



dg^- 



dN" 



K' 






Similarly, when we include the effects of correlations to order e in (III.ll) for a 
equation for N°- = F° becomes 

N'^it + At) = ^% {N\t) + e [IC\ - S\] N^{t) + IC''-^r>{t)} . 



(in.12) 
(in.13) 

{a}, the dynamical 
(IIL14) 



Note that if we set F" = in this equation, we get back the linearized Boltmann equation (IILS), since 






Inserting (III. 12), we can now write the exact equation of motion for the quantities A^" in the form of an 
infinite series. 



N^it + At) = A\ {N\t) + ej''^N^{t)) 



where 



J' 



Mb) 



x(c) 



J 



t(b) 



t(c) 



iC^iA^^IC^ 

a p 



A"^JC"JA''JC'' 

/9 7 '5 



■))• 



(III.15) 



(III.16) 



We now have an expression for the mean occupation number of a certain bit of the system at position x 
and time t + At, written as an infinite sum of terms, each of which is a function of the quantities TVg and Ni 
at nearby lattice sites x' and at previous time steps t' . As we consider terms in this series with more and 
more factors of AK,, the positions and times at which these quantities are evaluated will differ from x and 
t by greater amounts. However, for any given term in the series, the means Nq'^ ' {t') and N^^^'^ 'It') can 
be replaced by Nq it) and N'^ {t), and the expression ( III.16|) will only change by a quantity of order 
e, since spatial derivatives are ordered as e, and temporal derivatives are ordered as e^. Such a modification 
for a finite number of terms does not change the behavior of the system in the hydrodynamic limit. In 
fact, it follows that whenever the sum of terms in ( III.16| ) converges on a scale which goes to zero in the 
hydrodynamic limit, we can drop all the spatial and temporal variations in the single-particle means. Thus, 
(III. 15) can be rewritten in precisely the form of Eq. (III.3), where J is replaced by the renormalized matrix 



JMx) - E ^"'n(..) 



(III.17) 



yeL 



with all /C's in J^ evaluated at the point x and the time t. 



It is important to note that the above argument breaks down when the sum (III. 16) is divergent. In this 
case, t he effe cts of iarge scale variations in the iV^^'^'^'s must be considered. In general, for lattice gases 



where ( [11.16 ) is divergent, one must be quite careful about the analysis. For certain lattice gases, however, 
particularly systems in which the conserved quantities are ordered, we are interested in expanding around 
an equilibrium which is spatially invariant (for example, the FHP-I lattice gas). In this case, the zero-order 
means can be replaced by their universal values; however, one must still treat the spatial variation of the 
first-order means carefully. 

Now that we have rewritten the exact dynamical equation in a form commensurate with the original form 
of the lattice Boltzmann equation, the renormalized transport coefficients for the theory can be related to 
the eigenvalues of the matrix J'- in the same way that the original (Boltzmann) transport coefficients were 
related to the eigenvalues of the matrix J*^. Thus, if we can compute the matrix J*^ exactly, we can also 
compute the exact renormalized transport coefficients. 



We will now show that expression (p^II.16) for J can be written in a diagrammatic notation, allowing us 



to perform a perturbative calculation of J by summing over "Feynman diagram" -like objects, where the 
contribution from each diagram is just the product of factors associated with its vertices. The key observation 
which allows this reduction is the following: 

Theorem 1 For fixed a and P, JC^a can be broken down into a product of contributions from distinct vertices. 
Explicitly, 

x6Lq 

where the correlation vertex coefficients (CVC's), V, are defined by 

^%. = E(-i)'"''^''' ( n ^o] ( n ^Vt. (ni.i9) 

with the sum taken over ^J. ^ a^ and J^ ^ /5x, J' ^ -B. 

This theorem can be proven by expressing df°'/dT^\ and dg°'/dN^\ as products of the N^ and plugging 
into (III. 13). For a full derivation sea3. Using this result, it is possible to express every term in J'^j, in 



diagrammatic form. A general nonzero term is given by 

with a* and /3' fixed (i.e., not summed over), and |a'|, |/3'| > 2, except for the endpoints where a'^ ~ {a} 
and P^ = {b}. 

We define a diagram T by an integer k{T), which we refer to as the length of the diagram T, and a function 
aT(T), where for each t E {0, . . . , k} , axiT) Q B. Geometrically, we associate each a G arlT") with a virtual 
particle (VP) moving from {x{a),T) to (x(a) -I- c*("\r -t- 1) on the lattice A^+i = L x {0, . . . ,k + 1}. We 
refer to ariT) as the set of outgoing VP's for the diagram T. 

It is natural to define a corresponding set of incoming VP's for r > by I3t{t) = {b : a{i{b),x{b) — d'^''^) £ 
axir — 1)}. We also define ariT) = \aT{T)\ to be the total number of outgoing VP's for each value of r. 
Finally, given a diagram T, we can define a weight function 

w{T) = \{ n v-(%^(,)^. 



xS-L l<T<k{T) 

by taking the product of V over all vertices. 



The term ([11.20) can now be represented by the diagram T with aT{T) = a^ , where for consistency a" is 



defined to be the unique set with A o — 1- When r ^ 0, A!^ ^ — 1, so (3t{'t) — /3'^ for all r. It follows 



that the contribution from (III. 20) is exactly given by W{T). Thus, we can rewrite expression (III. 16) for 
J^";, as a sum over diagrams 



fc=iTer° (fe) 

where in general we define the set of diagrams T'^p^k) by 

r"^(fc) = {T : fc = fc(T), aril) > 1, for 1 < / < fc, arik) = a,/3T(l) = /?}. 

Note that V = (5 , so any diagram with incoming VP's at (x, t) but no outgoing VP's has weight zero. 
For many lattice gases certain other vertex factors V^^^ vanish also; diagrams with such vertices can be 



dropped from the sum (111.21). From (III. 17), J can now be written as a sum over diagrams in the same 
fashion as J^, 



k=lTeT\(x,k) 

where the set of diagrams to be summed over is given by 

r}(x,fc)^ U r'^(-i(fc). 



(III.22) 



IV. APPROXIMATIONS 



We have so far managed to write the exact formula for the hydrodynamic equations in the scaling limit 
only in terms of an infinite formal series. The natural next question to confront is whether this series can 
be summed. We would like to know whether the series is finite, and if we cannot sum the full series, at 
least we would like to find a set of reasonable approximations which we can make to truncate the series to 
one which is summable. The questions of convergence are rather difficult, and we will not address them 
here; in general, the convergence properties of the series depend on the form of the conserved quantities in 
the system. A variety of methods for performing partial sums of infinite series of diagrams while retaining 
physically importar± terms have been applied to related problems in quantum field theory and quantum 
many-body theoryBO. We will describe here several particular approximation methods which are useful for 
the kind of series which arise for known lattice gases. 



A. Short-r and SmalI-£ Truncations 



The simplest useful approximations involve truncating the sum (III. 21) to a finite number of terms by 
putting an upper bound on either the number of timesteps or the number of distinct nontrivial vertices 
allowed in each diagram. In the first case, the expression for the renormalized J-matrix is 



j(^);.(x) = ^ E w(n 

k=lTeT\{x,k) 



where the diagrams summed over are the same as those summed in Eq. (III. 22). Since for each fixed value of 
k there are a finite number of allowed diagrams, this sum is finite. We refer to this approximation as the short 
T approximation. In the second case, we allow k to be arbitrary, but allow only diagrams where the total 



number of vertices (x, k') with nonempty outgoing sets aT{k')y^ is less than or equal to some fixed number 
(.. We denote the sum restricted to these diagrams by J .. Again, there are only a finite number of such 
diagrams in this sum, which means that this sum must also be finite. This approximation is analogous to 
the weak-coupling expansions in quantum field theory, although in this case the coupling constants V^^ are 
usually not particularly small. The short-r and small-^ truncations give good consecutive approximations for 
many lattice gases, and in most instances pick out the contributions to the renormalization of the transport 
coefficients in decreasing order. In either of these two approximations, the Boltzmann approximation can be 
recovered, by taking r = I or £ = I . 

B. BBGKY Truncations 

Another good class of approximations, in which a reduced but still infinite set of diagrams is summed, 
corresponds to truncations of the BBGKY hierarchy of equations. Such an approximation involves neglecting 
q-CCF's with q > n for some fixed value of n. In our diagrammatic formalism, this amounts to restricting 
the sum to diagrams with axik') < n for 1 < k' < k. Whereas the total number of distinct diagrams in 
the complete sum grows exponentially in fc, the number of distinct diagrams in the BBGKY approximations 
grow polynomially, and are computationally more tractable. 

C. The Ring Approximation 

The n = 2 version of the BBGKY approximation is closely related to the ring approximation. The ring 
approximation is made by neglecting interactions between two propagating correlated quantities except at 
the initial and final vertices of a diagram. It is generally possible to calculate a closed-form expression for 
the infinite sum of diagrams corresponding to this approximation. Furthermore, it is usually fairly easy to 
calculate the asymptotic form of this approximation as fc — > cxd. This calculation often captures the most 
significant part of the long-term renormalization effects. In particular, for certain lattice gases which model 
two dimensional fluid systems, the ring approximation diverges logarithmically in \L\, which is in agreement 
with predictions from other theoretical frameworkaJ, and also with observed behaviouiil3. 

V. THE 1D3P LATTICE GAS 

As a simple example, we consider a diffusive lattice gas model in one dimension. The model has three 
bits per site (n = 3), corresponding to the presence or absence of left-moving, stationary, and right-moving 
particles, respectively. These bits are denoted by the respective elements of the set B = { — , 0, -I-}. Collisions 
occur only if exactly two particles enter a site. If we denote the two-particle states by + = { — , 0}, = { — , +}, 
and — = {0, +}, then the nontrivial elements of the state transition table can be written 



a 


{s -> s') 


s' 




+ 





- 


s 


+ 


1-nP 


nP{l - n'') 


nPjf 





nPn"- 


1-nP 


nP{\ - n"-) 


— 


nPil-n"-) 


nPn'- 


1-nP 



The bits nP and vT are stochastic bitsn which are sampled separately at each lattice site and at each timestep 
with average values (n^) = 2p and {pJ") = 1/2. Here, the parameter p G [0, 5] may be thought of as the 



^Note that r and p axe not indices here, but rather simply labels for the stochastic bits. 



probability of collision from, e.g., + to 0. The value of the bit n^ effectively determines whether or not a 
collision will occur, and that of n^ determines which of the two possible outcomes will result. 

Note that these collisions obey semi-detailed balance, since the columns of the above table sum to unity. 
Also note that they conserve particles {ric = 1); it is the particle density / that will obey the macroscopic 
diffusion equation. 

Beginni ng w ith t he ab ove state transition table, we can calculate the mean vertex coefficients V^, using 
equation ([II. 5) or (111.7). The nonzero mean vertex coefficients are given by 



vl = p + S\{l-3p). 

Note that since the ensemble-averaged collision operator is invariant under permutations (relabeling) on the 
bits, the mean vertex coefficients also have this symmetry. 

The nonvanishing correlation vertex factors are depicted graphically in Figure ll|; only a single vertex is 
shown in each equivalence class under the permutation symmetry. Note that the CVC's are also symmetric 
under an arbitrary permutation on the particle labels. 



l-2p/ 



-Pf 




V+ =-p/(l-/) 





VO =2p/(l-/) 





V°^ = -p 




V~^ = 2p 



V+^ = p(l - /) 



V-^ = 1 - 2p(l - /) 




V^„ = l 



FIG. 1. Vertex Factors for 1D3P Lattice Gas 
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The Chapman-Enskog analysis for this lattice gas is straightforward, and yields the following expression 
for the diffusivity in terms of the kinetic eigenvalue, A^: 



D^ 



-1 



3At V-A2 
The Boltzmann value for A^ is Spf. The renormalized result can be written 



(V.l) 



/ 



-3pf + 9p2/(l _ /) 



E 



A 
V 



E 



A 
M 



(V.2) 



V / 

where the notation in brackets indicates summation of the products of all internal ver tex factors over all 
diagrams with the depicted initial and final configurations. Together, Eqs. (V.l) and (V.2), with vertices 
given in Figure |l|, constitute an exact expression for the diffusivity of the 1D3P lattice gas. 

We have used a computer to numerically calculate the limit of the full fc-particle BBGKY approximation 
for certain values of / and p. The algorithm we used was to sum all diagrams of length < r on a lattice of 
width Z, then to take the limits as r, Z ^ oo. In this way, we have numerically approximated the limits of the 
/c-particle BBGKY sums for A: < 6. The results of this calculation, plotted as the correction to the Boltzmann 
value of the kinetic eigenvalue, versus mean occupation number, are graphed in Figure for p — 1/2. It 
is especially noteworthy that when / = 1/2 the ring approximation and lower BBGKY truncations are 
inadequate (indeed, the fc = 2 and ring truncations actually get the wrong sign for the correction), but the 
A: = 5 BBGKY truncation does quite well indeed. 



5X^ 




FIG. 2. Partial BBGKY Approximations for p = 1/2 



VI. CONCLUSIONS 



In this paper, we have outlined a complete kinetic theory of lattice gases, applied it to a simple model 
lattice gas, and compared the predictions of the theory to experiment. 
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The approach presented in this paper opens up a wide range of possible new work on discrete kinetic 
theory. By applying these techniques to compute deviations from the Boltzmann predictions for commonly 
used lattice gases, the results of simulations can be more accurately interpreted. Lattice gases are curpptly 
being used, both in industrial and academic settings, for computational fluid dynamics calculations^; to 
ensure the accuracy of these calculations, it is essential to account for the renormalization effects that we 
have studied here. 

In addition to quantitative refinement of lattice gas calculations, the theory presented here provides a tool 
with which to investigate fundamental physical phenomena in nonequilibrium statistical systems. In recent 
years, for example, lattice gases have been used to model many different hydrodynamic systems, includiae 
reaction-diffusion equations, and other systems capable of spontaneous self-organization. It is knowntHi 
that the Boltzmann approximation does not yield accurate results for the transport coefficients of such 
systems, unless the reactants are allowed to diffuse for several steps between reactions in order to artificially 
supress the correlations that developll3. Thus, it is likely that the methods described in this paper will be 
directly applicable to these systems, providing an important correction to their theory. More interestingly, 
these methods will also provide insight into the extremely subtle flow and agglomeration of interparticle 
correlations in pattern- forming lattice gases, and hence into the dynamical basis of self-organization. 
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